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ABSTRACT 

We have resolved the ultraluminous infrared galaxy (ULIRG), IRAS 09111 — 1007, 
with the new 350 /im-optimised Second Generation Submillimeter High Angular Res- 
olution Camera (SHARC II) and present the first submillimetre fluxes and images for 
the system. IRAS 09111 — 1007 comprises two interacting luminous infrared galax- 
ies (LIRGs) with a projected nuclear separation of 39 lij^ kpc. The Western galaxy is 
roughly four times more luminous in the submillimetre than its Eastern counterpart. It 
is an extremely bright LIRG with an AGN. The classification of the Eastern source is 
uncertain: it could be a Seyfert 2 galaxy or a LINER. We highhght IRAS 09111 - 1007 
as a system that necessitates further study: a double AGN ULIRG whose molecular 
gas content differs from other widely separated pairs and whose ULIRG phase might 
not be explained by current multiple merger and/or final stage ULIRG scenarios. 
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1 INTRODUCTION 

Amongst the first results of extragalactic mid-IR astronomy 
was the discovery of a small number of galaxies that emit the 
bulk of their bolometric luminosity in the infrared (Low & 
Kleinmann 1968, Kleinmann & Low 1970a,b). The InfraRed 
Astronomical Satellite, IRAS, detected large numbers of 
these ultraluminous infrared galaxies (ULIRGs) (Soifer et 
al. 1984, Joseph & Wright 1985, and Soifer, Neugebauer 
& Houck 1987) with quasar-like luminosities of Lir (8- 
1000 /im) > lO^'^ L©. There is still debate as to the nature 
of the far-IR power source in these galaxies: is the immense 
thermal energy driven by a dominant starburst, a dominant 
AGN or some combination of the two? These low-redshift 
IRAS-selected ULIRGs are expected to be the counterparts 
to the high redshift (2 > 1) SCUBA sources (see, e.g., Smail 
et al. 1998, Blain et al. 2002, Webb et al. 2003, Chapman et 
al. 2003). 

Most ULIRG systems have been shown to be disturbed, 
interacting or merging in some way when the separation 
of nuclei is less than lOkpc (Sanders et al. 1988, Clements 
et al. 1996, Murphy et al. 1996, Farrah et al. 2001). The 
nature of widely separated ULIRG systems is less clear 
(Dinh-V-Trung et al. 2001, Meusinger et al. 2001): are the 
components of the ULIRG (supposedly the end phase of 
the galactic interaction) beginning another merger or is the 



ULIRG a result of a multiple merger event (Borne et al. 
2000)? This latter scenario is possible in widely separated 
ULIRGs with resolved double nuclei but might not apply to 
the ULIRG system IRAS 09111 - 1007, which consists of 
two widely spaced but interacting luminous infrared galax- 
ies (LIRGs), each with a single nucleus. The two LIRGs 
have a projected separation of 39 h^^^ kpc and a velocity 
difference of 425 kms"^ (Due, Mirabel & Maza 1997). In 
this letter we present 350 /im resolved images and fluxes for 
IRAS 09111 - 1007. We model the far-IR dust emission to 
constrain the nature of the interaction. 



2 OBSERVATIONS AND DATA REDUCTION 

The data were taken using the Second Generation Sub- 
millimeter High Angular Resolution Camera (SHARC II) 
at the Caltech Submillimeter Observatory on Mauna Kea, 
Hawai'i, in January and March 2004. SHARC II is a 350 ^m- 
optimised camera (Dowell et al. 2003) built around a 12 x 32 
element close-packed bolometer array (Moseley et al. 2004). 
It achieves a point-source sensitivity of ~ 1 Jy Hz"^'''^ in 
good weather. The 384 pixels of the SHARC II array im- 
age a region of around 1.0' x 2.5'. Its filled absorber array 
provides instantaneous imaging of the entire field of view, 
sampled at roughly 2.5 pixels per nominal beam area. The 
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Name Coordinates 350/xm Flux log(LpiR.gED) log(Lpip|.ip(As) log(LiR) SFR 

{J2000) [Jy] [Lq] [Lq] [L0] [Mgyr-l] 

IRAS 09111 - 1007 11.91 11.86 12.09 81 



IRAS 09111 


- 1007W 


09 13 36.4 -10 19 31.8 


0.85±0.13 


11.80 


11.75 


11.98 


63 


IRAS 09111 


- 1007E 


09 13 38.8 -10 19 21.5 


0.23±0.04 


11.25 


11.20 


11.43 


18 



Table 1. Coordinates, 350/im fluxes, luminosities and star formation rates for the system and both components of IRAS 09111 — 1007 
(IRAS fluxes from Surace, Sanders &; Mazzarella 2004). The FIR (40 — 500 fiicn) and IR (8 — 1000 fim) luminosities are derived as follows; 
LpiRiSBD is computed from best-fitting single temperature SED, Lpjj^.jj^AS are calculated using the standard relations from 

Fullmer & Lonsdale (1989) and Sanders & Mirabel (1996) respectively. Note that the pair positions given in both NED and Surace, 
Sanders & Mazzarella (2004) are swapped. 



beam profile was measured on known compact sources, and 
was verified to be within five per cent of the difi'raction- 
limited beamwidth of 8.5". 

For these data the in-band zenith atmospheric opacity 
{tsso ij-m) ranged from 1.1 to 1.3, corresponding to a zenith 
transmission of around 30 per cent. Our observations were 
centred on the Eastern source of the pair, at position RA= 
09''13™'38.''6, Dec= -10°19'20" (J2000). 

In the submillimetre, emission from the atmosphere 
dominates the signal read out by the bolometers. In order 
to detect the faint celestial sources it is necessary to remove 
this emission. The sky signal, however, is largely correlated 
between pixels. At any moment in time, every pair of pixels 
can be used as a 'signal' beam and a 'reference' beam. Ad- 
ditionally, the telescope is moved such that each position on 
the sky can be viewed by many detectors at different times. 
This introduces a self-consistent self-calibration of all pixels 
compared with all others. A Lissajous scan pattern is used 
to ensure that the area on the sky is well covered and has 
substantial redundancy of observing. 

The data were reduced using the standard CSO reduc- 
tion software, CRUSH (Kovacs, in preparation). This soft- 
ware implements a self-consistent least-squares algorithm to 
solve for the celestial emission, taking into account instru- 
mental and atmospheric contributions to the signal. All ob- 
servations were taken using the Dish Surface Optimisation 
System (Leong 2003) , which corrects for the primary mirror 
deformation as a function of zenith angle, to improve the 
telescope efficiency and the pointing. 

The skymap is calibrated with a point spread function 
based on all point source observations (Callisto, Ceres & 
Arp 220) through the observing period at similar elevations. 
An oversampled fit is used to determine the position of 
the source and the flux per beam. The processed image is 
shown in Figure 1. 
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Figure 1. 350 fim continuum emission map of IRAS 09111 — 1007 
taken with SHARC II. The contours are in levels of 50 mjy/beam 
(dashed contours are negative and are an artifact of data reduc- 
tion). 



letter. The Eastern LIRG is called the 'Eastern source' or 
IRAS 09111 - 1007E (Murphy et al. 1996). Together they 
form the 'ULIRG system'. Figure 2 shows the resolved com- 
ponents of IRAS 09111-1007 imaged with the DSS, 2MASS 
and IRAS (HIRES processed; Surace, Sanders & Mazzarella 
2004) respectively. The 350 ^m fluxes for each component 
are presented in Table 1. The signal-to-noise in the detec- 
tion is 89 for the Western source and 29 for the Eastern 
source. The Western source is roughly four times more lu- 
minous in the submillimetre than its Eastern counterpart. 



3 RESULTS 

We are able to resolve the IRAS 09111 — 1007 system and ob- 
tain 350 nm fluxes and positions. An absolute pointing offset 
of 2.0" (Aa = 1.0", AS = 1.7") from the 2MASS positions 
was attributed to a pointing error (due to the imprecise ab- 
solute pointing knowledge of the CSO) and removed from 
the SHARC II image presented in this letter. 

The brighter LIRG, that which contributes «79 per cent 
of the total system 350 ^m emission, is known as the 'West- 
ern source' or IRAS 09111 — 1007W for the purposes of this 



4 SED MODELLING 

4.1 Dust Temperature Blackbody Fitting 

A modifled blackbody is used to model the total dust emis- 
sion of the IRAS 09111 - 1007 system, specifically: 

= {Mdust/D^)K{X)B,{\,Tdust) (1) 

where is the Planck function, k the mass absorption co- 
efficient of the dust («:(A) oc \~'^ , see Dwek 2004), Tdust and 
Mdust the equilibrium dust temperature and mass respec- 
tively, and D the distance of the galaxy. For the distance we 



IRAS 09111-1007 3 



I I I I I I I I I I I I I I I 



9 « 





113 39 38 




113 39 38 




I 13 39 38 



Figure 2. DSS Optical, 2MASS Ks, and IRAS 60 ^tm images of IRAS 09111 - 1007, with 350 ^tm contours overlaid on each (in levels of 
100 rnjy/beam). 
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Figure 3. Fit to the spatial variation in intensity for the HIRES- 
processed IRAS 60 lira image, cut through both galaxies. 



Model 


Source 




Tdust 

[K] 


Mdust 

[IO^Mq] 


Single Temperature 


System 
Western source 
Eastern source 


1.9 ±0.1 
1.9±0.1 
1.9 ±0.1 


31 ± 1 
31 ± 1 
31 ± 1 


220130 
170130 
5011° 


Two Temperature 

Mcold/M„arm = l 


System 


2.0 

(fixed) 


32 ± 2 
26 ± 5 


30011?° 


Two Temperature 

Mcold/Mwarm = ll 


System 


2.0 

(fixed) 


36 ± 3 
29 ± 1 


2801™ 


Single Temperature 
DL Dust Model 


System 
Western source 
Eastern source 


n/a 
n/a 
n/a 


30 ± 1 
30 ± 1 
30 ± 1 


25013° 
200130 
60+1° 



Table 2. Best-fitting SED parameters: /3, temperature and dust 
mass for the IRAS 09111 — 1007 system and components, with 
corresponding estimated Icr uncertainties. 



use WMAP cosmology (Hq = Tlkms^^ Mpc~\ = 0.27, 
and r^A = 0.73 (Bennett et al. 2003)). 

In order to better constrain the SED of each component, 
we determined the ratio of the fluxes of the two components 
at 60 /im. Beginning with the IRAS 60 /im HIRES-processed 
image of Surace et al. (2004) , we sliced along the axis of the 
known 350 /im sources. A pair of Gaussian intensity func- 
tions were fltted to the measured flux along this slice, con- 
straining them to have the same width and a flxed spacing 
as determined at 350 (see Figure 3) . The remaining four 
free parameters (one position, the width, and the two inten- 
sities) are then well-constrained, with the flux ratio being 
4.9 ± 0.8 (the total system flux as measured by IRAS is still 
valid, since the IRAS fluxes are derived with an aperture 
large compared to the source separation). The ratio derived 
from the 350 /xm measurement is 3. 7 ±1.1. Alternatively, the 
fraction of the flux from the brighter LIRG is 83 ± 3 per cent 
at 60 /^m and 79 ± 16 per cent at 350 /im, an entirely con- 
sistent measurement. The ratio of 79 per cent is accurate 
enough for the flts that follow. 

minimization determined the best-fitting emissivity 
index (/3) and corresponding single fit dust temperature for 
the system using the SHARC II 350 /im flux with the 60 /im 



and 100 /im IRAS fluxes derived in Surace et al. (2004) (the 
100 /im flux was used as the normalization). The best-fitting 
was /3=1.9±0.1, T=31±1K, and is shown in Figure 4. We 
also model both the Western and Eastern sources assuming 
that the IRAS fluxes are distributed in the same way as 
the 350 pim emission. Dust masses are calculated assuming 
the 100 ^m mass absorption coefficient, kioo, is 4Q era? 
(Draine & Lee 1984). 

With only two independent colours we cannot constrain 
the cold and warm dust masses and simultaneously fit the 
dust temperatures in a two temperature model. Using 13=2 
(as in Dunne & Eales 2001), if we assume equal amounts 
of cold and warm dust in the system we get best-fitting 
temperatures of Tcoid = 26 ± 5K and Twarm = 32 ± 2K. 
If we adopt a dust mass distribution typical of the highest 
luminosity galaxies from the SCUBA Local Universe Galaxy 
Survey (SLUGS), as given in Dunne & Eales (2001), we get 
a ratio Mcoid/Mwarm=ll. The best-fitting values are then 

Tcold = 29 ± 1 K, Twarm = 36 ± 3 K. 

As an alternative to these approaches, we also used a 
more physically-based dust grain model (the single temper- 
ature DL dust model - Draine & Lee 1984, Laor & Draine 
1993), an equal mixture of silicates and graphites with kioo 
of 31 and 54 cm^ g"^ for the silicate and graphite grains re- 
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IRAS 09111-1007 Best-fitting SEDs for /S=1.9, T=31K 
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Figure 4. Single Temperature Best-Fitting SED for the West- 
ern source (dashed), Eastern source (dot-dashed), and the 
IRAS 09111 - 1007 system (solid: (3=1.9, T=31K). 



spectively. In this case, the best-fitting temperature of the 
system was 30 ± 1 K. 

The dust temperature of ~ 31 K is within 1.5a of the 
average temperature of 38 ± 6 K from the previous 350 fim 
LIRG study by Benford (1999), with the first SHARC cam- 
era (although that survey did not include any widely sep- 
arated LIRG systems). The result that the 60 /xm dust is 
distributed in very similar ratios to the 350 ^im dust (Figure 
3) is responsible for the identical SED temperatures for both 
components. 

4.2 Ratio of Molecular Gas To Dust 

The molecular gas mass in the Western source is 2.3 x 
10^° M© (Mirabel et al. 1990). The Western source's dust 
meiss is derived from the best-fitting SED to give a molecu- 
lar gcis-to-dust ratio of 140 and 120 for the single temper- 
ature and DL dust models respectively. In the absence of 
21 cm Hi observations this value is a lower limit on the total 
(molecular -I- atomic) gas-to-dust ratio. 

The 8-1000 ^m luminosity (Lir, see Table 1) gives a 
value for the Western source IR luminosity-to-H2 mass ratio 
(LiR,/M(H2) - the star formation efficiency) of 42 Lq Mg^. 
The molecular gas-to-dust and IR luminosity-to-H2 mass 
ratios are consistent with values for the highest luminosity 
SLUGS galaxies (Dunne et al. 2000). 

Dinh-V-Trung et al. (2001) studied six widely separated 
(> 20 kpc) ULIRG systems in the complete 1 Jy sample of 
Kim & Sanders (1998). In their sample, the molecular gas 
was concentrated in the dominant source of the far-IR emis- 
sion. Although the Western source of IRAS 09111 - 1007 
is gas-rich, the Eastern source, by virtue of the amount of 
dust, is unlikely to be gas-poor. This would make this system 
different from their wide-pair ULIRG sample. 

4.3 Star Formation Rate 

The system and component luminosities (see Table 1) show 
that the two components form a ULIRG system, though 
neither is a ULIRG individually. The star formation rate 



in Table 1 uses the SED derived FIR luminosities and the 
relation of Thronson & Telesco (1986); 

SFR~*10-"' (^^)MQyr-' (2) 
assuming 'I' = 1 (typical values of ^' are 0.8-2.1). 



5 DISCUSSION 

5.1 Source Characterisation 

The optical image from the Digitized Sky Survey (Figure 2, 
left) shows the disturbed morphology of the Eastern source. 
In the 2MASS Ks band image (Figure 2, centre) the West- 
ern source is slightly more luminous than the Eastern, a 
trait even more prominent in the submillimetre (Figure 1). 
Merger models such as Barnes & Hernquist (1991) predict 
gas and dust to be concentrated in the galaxy centre as the 
ULIRG interaction condenses large amounts of the ISM into 
the nuclear region. 

Line ratios from Due et al. (1997) classify the Western 
source as a Seyfert 2. Dudley (1999) found prominent poly- 
cyclic aromatic hydrocarbon features in the 8-13 fim dust 
emission spectra - indicative of a starburst. The Eastern 
source is either a Seyfert 2 or LINER galaxy (Due et al. 
1997), while Gongalves, Veron-Cetty & Veron (1999) also 
find evidence for a starburst. Neither of the two sources in 
the IRAS 09111 - 1007 system were detected in the ROSAT 
All-Sky Survey. A non-detection in the ROSAT band does 
not necessarily mean that a source is intrinsically weak since 
the soft X-ray band is sensitive to X-ray absorption, which 
is common in AGN. 

5.2 Merging Stage 

Without additional submillimetre data we are unable to con- 
strain the relative temperatures and masses of the cold and 
warm dust components. Whether the dust temperature of 
the system is related to the stage of merging is not clear. 
Mazzarella, Bothun & Boroson (1991) find an increase in 
warm dust temperature with merging stage, although Klaas 
et al. (2001) argue the cold dust temperature would increase 
as well. The wide separation of the pair would suggest that 
IRAS 09111 — 1007 is at the beginning of a merger, a notion 
supported by the value of the IR luminosity-to-H2 mass ra- 
tio, which falls within a region of Lir/M(II2) vs Lir space 
that is common for early merging systems (Sanders, Scoville 

6 Soifer 1991). 

5.3 Widely Spaced ULIRG Pairs 

Unlike the wide pair sample of Xu & Sulentic (1991) both 
components are enhanced in the far-IR. With a velocity dif- 
ference of 425kms~^ (Due et al. 1997) and a projected sepa- 
ration of 39 h^^ kpc, Monte Carlo simulations give the prob- 
ability of the pair being bound as 0.88 (Schweizer 1987). 

Although widely separated ULIRG pairs are not uncom- 
mon the nature of their interaction is still uncertain. Murphy 
et al. (1996) postulated the presence of a third nucleus in 
widely separated pairs, though no double nucleus has been 
detected in either component of IRAS 09111 - 1007. The 
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galaxies were shown to be unresolved at 0."5 resolution in 
a 6 cm search by Crawford et al. (1996). However the non- 
detection of a double nucleus cannot rule out a multiple 
merger since the time-scale of nuclei coalescence is short 
(Surace, Sanders & Evans 2000, Meusinger et al. 2001). In 
multiple approach merger models (e.g., Dubinski, Mihos & 
Hernquist 1999) the merging process is a series of encoun- 
ters where bound components approach and separate. A pre- 
vious encounter may have triggered the starburst/AGN in 
the system. High resolution optical imaging (Borne et al. 
2000) could decide between these scenarios by either detect- 
ing multiple nuclei or confirming single nuclei. High reso- 
lution CO imaging would be needed to detect whether gas 
has been disturbed by a previous phase of the merging event 
(Mihos & Hernquist 1996). 



6 CONCLUSIONS 

We have resolved the widely separated ULIRG system of 
IRAS 09111 - 1007 with the SHARC II detector at 350 /im. 
This system comprises two LIRGs with a projected sepa- 
ration of 39 kpc, or around two optical diameters. The 
Western component dominates the far-IR flux at both 60 
and 350 /xm, carrying 79 per cent of the total system lumi- 
nosity. Although the luminosity of the system is large, our 
fluxes suggest a dust temperature of 31 K for this system, 
with both components at the same temperature to within 
the sensitivity of this measurement. The wide separation 
and the value of the IR luminosity-to-H2 mass ratio sug- 
gest that the pair are at an early stage of interaction. But 
the high luminosity of the system (Lir = 1.2 x IO^'^Lq) 
would be unusual for such a stage, unless the components 
had experienced a previous merger or interaction. A high 
resolution optical search for multiple nuclei within each com- 
ponent is needed. Their absence could indicate that the dou- 
ble AGN-LIRG system of IRAS 09111 - 1007, and perhaps 
other widely spaced ULIRG pairs, might be unexplained by 
current theories of ULIRG formation and evolution. 
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